Absence of a day-night asymmetry in the ^Be solar neutrino rate in Borexino 
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Abstract 

We report the result of a search for a day-night asymmetry in the ^Be solar neutrino interaction rate in the Borexino 
detector at the Laboratori Nazionali del Gran Sasso (LNGS) in Italy. The measured asymmetry is A^,, - 0.001 ± 0.012 
(stat) ± 0.007 (syst), in agreement with the prediction of MSW-LMA solution for neutrino oscillations. This result 
disfavors MSW oscillations with mixing parameters in the LOW region at more than 8.5 cr. This region is, for the first 
time, strongly disfavored without the use of reactor anti-neutrino data and therefore the assumption of CPT symmetry. 
The result can also be used to constrain some neutrino oscillation scenarios involving new physics. 
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In the last two decades solar neutrino |[Tl|2l[3l and re- ? 

actor anti-neutrino |4| experiments have demonstrated a 

that solar electron neutrinos undergo flavor conversion g 

along their trip from the Sun's core to the Earth. The io 

conversion is well described by the so-called Mikheyev- n 

Smirnov-Wolfenstein (MSW) matter-enhanced neutrino 12 



oscillations 15] with Large Mixing Angle (LMA) oscil- 
lation parameters. A generic feature of matter-enhanced 
neutrino oscillations is the potential for the coherent re- 
generation of the flavor eigenstate when solar neutri- 
nos propagate through the Earth |6|, as they do during 
the night. Thus, there is the potential for those solar neu- 
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trino experiments that are principally, or entirely, sensi- 
tive to Ve to detect different solar neutrino interaction 
rates during the day and during the night. Solar neu- 
trino day-night asymmetry measurements are sensitive 
to both Vf appearance and disappearance. 

The magnitude of this day-night effect is expected to 
depend on both neutrino energy and the neutrino os- 
cillation parameters. Previous experiments |7, 8| have 
shown that for high energy (~5-15 MeV) solar neutri- 
nos, the day-night asymmetry is less than a few per- 
cent, in agreement with the MSW-LMA prediction. 
At lower neutrino energies (around 1 MeV), the pre- 
dicted day-night asymmetry for MSW-LMA is also 
small (<0.1%) [9 1; however, other scenarios including 
different MSW solutions and neutrino mixing involv- 
ing new physics |10| predict much larger day-night 
effects. For example, in the so-called LOW region 
(10^^ eV^<6tn'^<lO-^ eV^) of MSW parameter space, 
which is currently strongly disfavored only by the Kam- 
LAND anti-neutrino measurement under the assump- 
tion of CPT symmetry, the day-night asymmetry would 
range between about 10% and 80% for neutrino ener- 
gies near 1 MeV. We present here the first measurements 
sensitive to the day-night asymmetry for solar neutrinos 
below 1 MeV. This result is an essentially new and inde- 
pendent way to probe the MSW-LMA prediction and is 
potentially sensitive to new physics affecting the propa- 
gation of low energy electron neutrino in matter. Partic- 
ularly, this result is independent from the KamLAND 
measurement, which probes anti-neutrino interactions 
at higher energies (>1.8 MeV). 

The Borexino experiment at LNGS detects low en- 
ergy solar neutrinos by means of their elastic scattering 
on electrons in a large volume hquid scintillator detec- 
tor Real-time detection (with fis absolute time reso- 
lution) of all events is made by collecting the scintilla- 
tion light with a large set of photomultipliers. The very 
low intrinsic radioactivity of the scintillator and of the 
materials surrounding it allows a clean spectral separa- 
tion between the neutrino signals and the residual back- 
ground. As the neutrino-electron elastic scattering cross 
section is different for and v^-v^, Borexino can mea- 
sure the electron neutrino survival probability and is, as 
a result, sensitive to the day-night effect. 

We recently released a precise measurement of the 
^Be neutrino interaction rate in Borexino with a total 
uncertainty less than 5% (TT]. In this Letter, we present 
a study of the day-night asymmetry in the same ^Be so- 
lar neutrino rate, placing a stringent limit on the size of 
the possible effect. We show that this limit improves 
the constraint on the solar neutrino oscillation parame- 
ters from solar neutrino experiments alone and excludes 
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Figure 1 : The experimental exposure function (black continuous line) 
and the ideal exposure function (red dotted line). The interval from 
-180° to -90° corresponds to day time and the one from -90° to 0° 
to night time. We recall that at LNGS latitude the Sun is never at the 
zenith. 



new physics scenarios that cannot be rejected with the 
currently available data. 

The Borexino detector lfT2l [131 is located in Hall C 
of the Laboratori Nazionali del Gran Sasso (latitude 
42.4275° N) in Italy and has taken data since May 2007. 
The sensitive detector consists of ~278 tons of very pure 
organic liquid scintillator contained in a 4.25 m radius 
nylon vessel. The scintillator is viewed by 2212 photo- 
multipliers and is shielded against external neutrons and 
-y radiation llT4ll . The energy of each candidate event is 
measured by the total amount of collected light, while 
the position of the event is reconstructed using the time- 
of-flight of the light to the photomultipliers. 

The data used in this analysis were collected between 
May 16"', 2007 and May 8"', 2010 and correspond to 
740.88 live days after applying the data selection cuts. 
We define "day" and "night" using 0-, the angle between 
the vertical z-axis of the detector (positive upward) and 
the vector pointing to the detector from the Sun, follow- 
ing 1 2]. Note that, with this definition, cos 6, is negative 
during the day and positive during the night. The dis- 
tance that the neutrinos propagate within the Earth is 
small for negative cos0, (the ~1.4 km LNGS overbur- 
den) and ranges up to 12049 km for positive cos 0-. Our 
day and night livetimes were 360.25 and 380.63 days, 
respectively. The distribution of 0~ corresponding to the 
live time (experimental exposure function) is shown in 
Fig. [T] and its asymmetry with respect to -90° is mainly 
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due to maintenance and calibration activities which are 
normally carried out during the day. 

As discussed in UJJ, scintillation events due to ^Be 
solar neutrinos cannot be distinguished from back- 
ground events (cosmogenics and radioactivity) on an 
event-by-event basis. The signal and background con- 
tributions are therefore determined using a spectral fit to 
the energy spectrum of the events reconstructed within a 
suitable fiducial volume (86.01 m^^ in [11 J ), and passing 
a series of cuts which eliminate muons and short-lived 
cosmogenic events, time correlated background events, 
and spurious noise events (details of this event selection 
will be published in ifTSll ). The experimental signature 
of the mono-energetic 862 keV ^Be solar neutrinos is a 
Compton-like electron scattering "shoulder" at approx- 
imately 660 keV. 

In the analysis reported here we use a spherical fidu- 
cial volume significantly larger than the one used in 
ifTTIl in order to increase the size of the data sample. This 
choice is justified by the fact that the additional exter- 
nal background that enters this larger fiducial volume is 
due to gamma radioactivity emitted by the materials sur- 
rounding the scintillator volume. As this background is 
expected to be the same during day and night, it should 
not affect the day-night asymmetr}{^ 

We determined that our sensitivity to the day- 
night effect is maximized by a 3.3 m fiducial ra- 
dius, which gives a 132.5 ton fiducial mass containing 
4.978x10^' e". With this choice of fiducial mass, the 
signal-to-background (S/B) ratio in the "^Be neutrino 
energy window" (550 to 800 keV) is 0.70 ± 0.04. This 
value is smaller than the one in fTTl due to the increase 
in spatially non-uniform backgrounds produced by ex- 
ternal gamma rays and ^^-Rn events. 

The day-night asymmetry, A^/,,, of the ^Be count rate 
is defined as: 



Rn + Rd ~ Jr) 



(1) 



where 7?^ and R^, are the ^Be neutrino interaction rates 
during the night and the day, respectively, Rdiff is their 
difference, and {R} is their mean. 

Fig. |2] shows the day and night energy spectra super- 
imposed and normalized to the same live-time (the day 
one), while Fig.[3]shows the 9- distribution of the events 
in the ^Be neutrino energy window normalized by the 
experimental exposure function. By using the total ^Be 



As explained later on, not all backgrounds relevant for this anal- 
ysis are the same during day and night. Particularly, the background 
induced by ^'"Po as is not the same because of the long ^'"Po lifetime 
and of the different length of days and nights in summer or winter. 
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Figure 2: The energy spectrum of events during day (red) and night 
(black) normalized to the day live-time in the enlarged FV. The insert 
shows the 'Be neutrino energy window. See 1111 for details on this 
spectral shape. 



count rate measured in [llj, a correction has been ap- 
plied to the exposure function to account for the annual 
modulation of the neutrino flux due to the seasonal vari- 
ation of the Earth-Sun distance. Before correction, the 
asymmetric distribution of our day and night livetime 
throughout the year is expected to increase the measured 
^Be neutrino count rate by 0.37% during the night and 
decrease it by 0.39% during the day. The day and night 
spectra in Fig |2]are statistically identical, as proved by 
the fit to the data shown in Fig. [3] Indeed, by fitting 
with a constant distribution the data in Fig.|3]we obtain 
a probability = 0.44. Any deviation from a straight 
line would be a signature of day-night modulation. For 
illustration, we include in Fig. [3] the expected shape for 
the LOW solution (Am^^^l.O ■ 10"^ eV^ and tan2(6'i2) 
- 0.955). Fitting the distribution with a flat straight 
line yields ;t'^/ndf = 141.1/139, showing that the data 
are consistent with the no day-night effect hypothesis. 

One way to quantitatively constrain A^„ is to deter- 
mine Ro and R^ separately by independently fitting 
the day and night spectra using the same spectral fit- 
ting technique used in determining the total ^Be flux 
in ifTTIl and then comparing the results using Eq. 1. 
Note that because these neutrinos are mono-energetic, 
we expect the shape of the ^Be electron recoil spec- 
trum to be identical during day and night. This yields 
Ad„ = 0.007 + 0.073. This method has the virtue of al- 
lowing for the possibility of different background rates 
during day and night. However, this analysis is less sen- 
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Figure 3 : Normalized 0; -angle distribution of the events in the FV in 
the 'Be neutrino energy window. The effect of the Earth's elliptical 
orbit has been removed. The blue line is the expected effect with the 
LOW solution (Am;, = 1.0 ■ 10"' eV^ and tan-(6>i2) = 0.955). 



sitive than the one described below and is not used for 
the final result. 

A stronger constraint on A^/„ is obtained by mak- 
ing the very reasonable assumption that the main back- 
grounds that limit the sensitivity in LI 1 1 (**^Kr and ^"^Bi) 
are the same during day and night. With this assump- 
tion, Arfn is obtained by subtracting the day and night 
spectra (normalized to the day live time) following the 
second term in Eq. ([TJ and then searching for a resid- 
ual component having the shape of the electron recoil 
spectrum due to ^Be neutrinos. If Arf„=0 and the back- 
ground count rates were constant in time the subtracted 
spectrum would be flat. 

The subtracted spectrum is shown in Fig.|4] where the 
lower plot is a zoom of the upper one in the energy re- 
gion between 0.55 and 0.8 MeV. The result is a flat spec- 
trum, consistent with zero, except for a clear negative 
^"'Po peak visible in the low energy region. This nega- 
tive peak arises because the ^"^Po background count rate 
in Borexino is decaying in time (ti/2 - 138.38 days), 
and the day and night livetime are not evenly distributed 
over the 3 years of data taking. The ^^"Po count rate was 
highest at the time of the initial filling in May 2007, and 
has since decayed. Therefore, the ^"'Po count rate has 
been higher on average during the summers (when days 
are longer), leading to a noticeable effect in the sub- 
tracted spectrum. This effect is taken into account by 
including both the ^'"Po and ^Be spectral shapes in the 
fit. Fitting between 0.25 and 0.8 MeV, we obtain R^ff - 
0.04+0.57 (stat) cpd/100 t. The amplitude of the result- 
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ing electron recoil spectrum induced by the interaction 
of ^Be neutrinos is too small to be shown in Fig. |4] In 
order to see its spectral shape we plot the recoil spec- 
trum with an amplitude corresponding to the expected 
day-night asymmetry for the LOW solution. 

The R^iff result is confirmed by removing alpha events 
from the day and night spectra using a pulse shape anal- 
ysis based statistical subtraction technique [11] before 
creating the difference spectrum. In this case, no resid- 
ual ^'"Po peak is expected or observed in the difference 
spectrum. Fitting the data between 0.25 and 0.8 MeV 
using only the ^Be recoil shape yields a result consis- 
tent with the previous one. The difference in the central 
values is included in the systematic uncertainty. 

Using (R) =46+1.5 (stat) +[-^ (syst) cpd/100 t CT] 
we obtain Ad„ = 0.001 ± 0.012 (stat) ± 0.007 (syst) from 
Eq. 1 . The statistical error in Ad„ is given by 



{R) \ 



(^diff (r\<R>)\ 



n2 



<R>^ 



o'(Rdiff) 
(R) 



because the total relative experimental error associated 
with (R) is negligible with respect to ""^''^^^^ . 

The main systematic errors are listed in Table [T] The 
dominant uncertainties are associated with the differ- 
ence between the R^ff central values obtained with and 
without statistical subtraction of the a events, and the 
maximum effect on TJjiff from potential small changes 
in the '"'Bi background in the detector. These uncer- 
tainties will be detailed in 1 15 |. 

This new tight constraint on the day-night effect in 
^Be solar neutrinos has interesting implications on our 
understanding of neutrino oscillations. To investigate 
this, we calculated the expected day-night asymme- 
try for 862 keV neutrinos under different combinations 
of mixing parameters in the MSW oscillation scenario. 
The comparison of these predictions with our experi- 
mental number is displayed on the right panel of Fig.|5] 
The red region is excluded at 99.73% c.l. (2 d.o.f.). In 
particular, the minimum day-night asymmetry expected 
in the LOW region (10^^ eV^ < Am^ < 10"^ eV^) is 



Source of error 


Error on A^„ 


Live-time 
Cut efficiencies 
Variation of ^'°Bi with time 
Fit procedure 


<5-10-* 
0.001 
+0.005 
± 0.005 


Total systematic error 


0.007 



Table 1 : List of systematic errors on Aj„ . 
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Figure 4: Difference of niglit and day spectra in tlie FV. Tiie fit is 
performed in the energy region between 0.25 and 0.8 MeV witii tlie 
residual ^'"Po spectrum and the electron recoil spectrum due to the 
'Be solar neutrino interaction. The fit results are in cpd/100 1. The top 
panel shows an extended energy range including the region dominated 
by the ' ' C background while the bottom panel is a zoom of the ■'Be 
energy window between 0.55 and 0.8 MeV. The blue curve shows the 
shape of electron recoil spectrum that would be seen assuming the 
LOW solution as in Fig. |3] 



0.117, which is more than 8.5 cr away from our mea- 251 

surement, assuming gaussian errors for A,/„. 252 

This effect can also be seen in a global analysis of all ^53 

solar neutrino data. We have carried out such an analy- 254 

sis, assuming two neutrino oscillations (i.e. fl^ = 0, we 255 

have checked that the inclusion of the third family does 256 

not change any of the conclusions and will be published 257 

in 1 1 5 1 ), including the radiochemical data [ 1 ] , the Super- 258 

Kamiokande phase 1 and phase III data ||2l, and the SNO 259 

LETA data and phase 111 rates f3\. The analysis takes 260 

into account the experimental errors (the systematic and 261 

statistical errors summed in quadrature) and the theoret- 262 

ical errors in the total count rates, including the correla- 263 

tion of the ^Be and '^B theoretical fluxes |16|. We use 264 

flux predictions from a recent high metallicity standard 265 

solar model flT^ and we include the bin-to-bin corre- 266 

lations in the uncertainties in the predicted ^B neutrino 267 

recoil spectrum resulting from the uncertainties in the 268 

predicted neutrino spectrum, and from energy threshold 269 

uncertainties and energy resolution in the experiments. 270 

The left panel of Fig. |5] shows the 68.27, 95.45 and 271 
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Figure 5: Neutrino oscillations parameter estimation in three solar 
neutrino data analyses (with 2 d.o.f.): 1) 99.73% c.l. excluded region 
by the Borexino 'Be day-night data (hatched red region in the right 
panel); 2) 68.27%, 95.45%, and 99.73% c.l. allowed regions by the 
solar neutrino data without Borexino data (left panel); 3) Same c.l. 
allowed regions by all solar neutrino data including Borexino (filled 
contours in right panel). The best fit point in the left (right) panel 
is Am^ = (5.2 ^l^p -IQ-K tan^e=OAl (0.46 ^jj-jj*). The LOW 
region is strongly excluded by the 'Be day-night data while the al- 
lowed LMA parameter region does not change significantly with the 
inclusion of the new data. 



99.73% c.l. neutrino mixing parameter regions allowed 
by all solar neutrino data without Borexino. The best-fit 
point is in the LMA region (Am^ - (5.2 -10"^ eV^ 
and tan^ 6 =0.47+° °^) and a small portion of the LOW 
region is still allowed at A^^ = 11. 83. 

The right panel of Fig. |5] shows the regions of al- 
lowed parameter space after adding the Borexino data 
(the ^Be total count rate ifTTl , the day-night asymme- 
try reported in this paper, and the **B total count rate 
above 3 MeV (0.22 + 0.04 (stat) ± 0.01 (syst)) cpd/100 
t and spectral shape (5 bins from 3 to 13 MeV) |18|) to 
the analysis. The LMA region is only slightly modified 
(the new best fit point is Am^ - (5.2 +Q g) 10"^ eV^ and 
tan^ e - 0.46+0 03), but the LOW region is strongly ex- 
cluded at Ax^ > 190. Therefore, after the inclusion of 
the Borexino day-night data, solar neutrino data alone 
can single out the LMA solution with very high confi- 
dence, without the inclusion of anti-neutrino data and 
therefore without invoking CPT symmetry. 

This result is an essentially new and independent way 
to probe the MSW-LMA prediction and is potentially 
sensitive to new physics affecting low energy electron 
neutrino interactions. As an example, we note that our 



272 day-night asymmetry measurement is very powerful in 

273 testing mass varying neutrino flavor conversion scenar- 

274 ios. We find, for example, that our A^n data excludes 

275 the set of MaVaN parameters chosen in IfTOl to fit all 

276 neutrino data at more than 10 cr. 

277 In conclusion, we have searched for a day-night 

278 asymmetry in the interaction rate of 862 keV ^Be so- 

279 lar neutrinos in Borexino. The result is A^,, = 0.001 + 

280 0.012 (stat) ± 0.007 (syst), consistent both with zero and 

281 with the prediction of the LMA-MSW neutrino oscilla- 

282 tion scenario. With this result, the LOW region of MS W 

283 parameter space is, for the first time, strongly disfavored 

284 by solar neutrino data alone. The result constrains cer- 

285 tain flavor change scenarios involving new physics. 
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